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Abstract 

Macroscopic quantum coherence and spin-phase interference are studied be- 
tween excited levels in single-domain ferromagnetic particles in a magnetic 
field along the hard anisotropy axis. The system has the general structure of 
magnetocrystalline anisotropy, such as one showing a biaxial, trigonal, tetrag- 
onal, and hexagonal symmetry. This study not only just yields the previous 
spin-phase interference results for the ground state tunneling, but also pro- 
vide a generalization of the Kramers degeneracy to coherently spin tunneling 
at low-lying excited states. These analytical results are found to be in good 
agreement with the numerical diagonalization. We also discuss the transition 
from quantum to classical behavior and the possible relevance to experiment. 
PACS number(s): 75.45. -Fj, 75.50.Jm, 03.65.Bz 







I. INTRODUCTION 



Macroscopic quantum phenomena (MQP) in nanoscale magnets have received much at- 
tention in recent years both from theories and from experiments-lll A number of nanoscale 
samples in different systems have been identified as the promising candidates for the obser- 
vation of macroscopic quantum tunneling (MQT) and coherence (MQC). Maybe even more 
interesting subject is that the topological Wess-Zumino term, or Berry phasell can lead to 
remarkable spin-parity effects for some spin systems with high symmetries.0iHi It has been 
shown that the ground-state tunnel sphtting is completely suppressed to zero for half-integer 
total spins in biaxial ferromagnetic (FM) particles in the absence of a magnetic field due to 
the destructive phase interference between topologically different tunneling paths.! However, 
the phase interference is constructive for integer spins, and hence the splitting is nonzero. 
The spin-parity effects can been interpreted as Kramers' degeneracy at zero magnetic field, 
but in the case of a field along the hard anisotropy axis, these effects are not related to 
the Kramers' theorem since the field breaks the time reversal symmetry.0i Experiment on 
Fes showed that the oscillation of the tunnel splitting as a function of the magnetic field 
along the hard anisotropy axis was due to quantum interference of two tunnel paths with 
opposite windings, which was a direct evidence of the topological part of the quantum spin 
phase (Berry phase) in a magnetic system.0 Recent theoretical and experimental studies 
include the quantum relaxation in magnetic molecules,! the spin tunneling in a swept mag- 
netic field,0 the thermally activated resonant tunneling with the help of the perturbation 
theory0 and the exact diagonalization,lli the auxiliary particle method,lii the discrete WKB 
method and a nonperturbation calculation,0 the non-adiabatic Landau-Zener model,lll the 
calculation based on exact spin-coordinate correspondence, 111 and the effects caused by the 
higher order term and the nuclear spins on the tunnel splitting of Fes.S 

Up to now theoretical studies have been focused on spin-phase interference at excited 
levels in simple biaxial FM particles! or at ground states in FM particles with general 
structure of magnetocrystalline anisotropy.0 However, the spin-phase interference between 
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excited-level tunneling paths is unknown for FM particles with a general structure of mag- 
netocrystalline anisotropy. The purpose of this paper is to extend the previous results to 
resonant quantum tunneling and spin-phase interference at excited levels for single-domain 
FM particles in the presence of a magnetic field along the hard anisotropy axis. Moreover, 
the system considered in this paper has a general structure of magnetocrystalline anisotropy, 
such as biaxial, trigonal, tetragonal and hexagonal symmetry around which has two, three, 
four and six degenerate easy directions in the basal plane at zero field. Therefore, our study 
provides a nontrivial generalization of the Kramers degeneracy to coherently spin tunneling 
at ground states as well as low-lying excited states in a magnetic field. 

To compute the tunnel splitting, we consider the imaginary time transition amplitude 
in the spin-coherent-state path-integral representation. Integrating out the moment in the 
path integral, the spin tunneling problem is mapped onto a particle moving problem in one- 
dimensional periodic potential U (0). By applying the periodic instanton method, we obtain 
the low-lying tunnel splittings of the nth degenerate excited states between neighboring po- 
tential well. The periodic potential U (0) can be regarded as a one- dimensional superlattice. 
The general translation symmetry results in the energy band structure. By using the Bloch 
theorem and the tight-binding approximation, we obtain the low-lying energy level spectrum 
of the excited states. Our results show that the tunnel splittings depend significantly on the 
parity of the total spins of FM particles. External magnetic field yields an additional contri- 
bution to the Berry phase, resulting in oscillating field dependence of the tunnel splittings 
for both the integer and half-integer total spins. These analytical results are found to be in 
good agreement with the exact diagonalization computation. And the structure of energy 
level spectrum for the trigonal, tetragonal and hexagonal symmetry is found to be much 
more complex than that for biaxial symmetry. The transition from quantum to classical 
behavior is also studied and the second-order phase transition is shown. Another important 
conclusion is that the spin-parity effects can be reflected in thermodynamic quantities of 
the low-lying tunneling levels. This may provide an experimental test for the spin-parity or 
topological phase interference effects in single-domain FM nanoparticles. 
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II. THE PHYSICAL MODEL 



For a spin tunneling problem, the tunnel splitting for MQC or the decay rate for MQT 
is determined by the imaginary-time transition amplitude from an initial state \i) to a final 
state I/) as 



where VQ, = sin OdOdcj) is the measure of the path integral. For FM particles at sufficiently 
low temperature, all the spins are locked together by the strong exchange interaction, and 
therefore only the orientation of magnetization M can change but not its absolute value. In 
the spin-coherent-state representation the Euclidean action Se can be written as 



where Mq = |M| = hjS/V, S is the total spins, V is the volume of the particle, and 7 is the 
gyromagnetic ratio. The first two terms in Eq. (2) define the Wess-Zumino term (or Berry 
phase) which arises from the nonorthogonality of spin coherent states. The Wess-Zumino 
term has a simple topological interpretation. For a closed path, this term equals —iS times 
the area swept out on the unit sphere between the path and the north pole. The first term 
in Eq. (2) is a total imaginary-time derivative, which has no effect on the classical equations 
of motion, but it is of crucial importance for the spin-parity effects. 

It is noted that from 6Se = of the action Eq. (2) reproduces the classical equation 
of motion whose solution is known as an instanton, and describes the (1 © l)-dimensional 
dynamics in the Hamiltonian formulation, which consists of the canonical coordinates and 
Pfj, = S{l — cos6). According to the standard instanton technique in the spin-coherent- 
state path-integral representation, the tunneling rate F for MQT or the tunnel splitting A 
for MQC is given byEI 




(1) 




cose + E{e,^) , 



(2) 




(3) 
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where ujp is the oscillation frequency in the well, Sd is the classical action, and the prefactor 
A originates from the fluctuations about the classical path. It is noted that Eq. (3) is based 
on quantum tunneling at the level of ground state, and the temperature dependence of the 
tunneling frequency (i.e., tunneling at excited levels) is not taken into account. The instanton 
technique is suitable only for the evaluation of the tunneling rate or the tunnel splitting at the 
vacuum level, since the usual (vacuum) instantons satisfy the vacuum boundary conditions. 
Recently, different types of pseudoparticle configurations (periodic or nonvacuum instantons) 
are found which satisfy periodic boundary conditions.! 

For a particle moving in a double- well-like potential U (x) , the WKB approximation gives 
the tunnel splitting of the nth excited levels a 

Aij„ = ^^^exp[-5(i?„)], (4) 

71 



and the imaginary-time action is 



S {En) = 2V2m / dxJU (x) - E„, (5) 

where Xi^2{En) are the turning points for the particle oscillating in the inverted potential 
—U (x). oj (En) = 2n/t (En) is the frequency of oscillations at the energy level En, and t (En) 
is the period of the real-time oscillation in the potential well, 

t {En) = V2^ / , 6 

■JxsiE^) JE„-U(x) 



where 0:3^4 {En) are the classical turning points for the particle oscillating inside U {x). The 
functional-integral and the WKB method show that for the potentials parabolic near the 
bottom the result Eq. (4) should be multiplied by ^ .BH This factor approaches 

1 with increasing n and it is very close to 1 for all n: 1.075 for n = 0, 1.028 for n = 1, 
1.017 for n = 2, etc. Stirling's formula for n\ shows that this factor trends to 1 as n ^ cxd. 
Therefore, this correction factor, however, does not change much in front of the exponentially 
small action term in Eq. (4). 
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III. MQC FOR BIAXIAL SYMMETRY 



In this section, we consider an FM system with biaxial symmetry in a magnetic field 
along the hard anisotropy axis. The magnetocrystalline anisotropy energy can be written as 

E {6, 0) = Ki cos^ e + K2 sin^ 6 sin^ - MqH cos 6 + Eq, (7) 

where Ki and K2 are the longitudinal and the transverse anisotropy coefficients satisfying 
Ki ^ K2 > 0. Eq is a constant which makes E{6, (p) zero at the initial orientation. Although 
the tunnel splittings and spin-phase interference effects of this system can be easily obtained 
by direct numerical diagonalization of the Hamiltonian (see Refs. 4 and 7, and in the 
following), it is of interest to understand these features analytically. 
Adding some constants, we rewrite Eq. (7) as 

E{6, (p) — Ki (cos 6 — cos sin ^sin 0, (8) 

where cos 6*0 = ^^j^- As Ki ^ K2 > , the magnetization vector is forced to lie in the — 9q 
plane, and therefore the ffuctuations of 9 about ^0 are small. Introducing 9 — 9Q-\-a {\a\ <g 1), 
the total energy E [9, 0) reduces to 

E (a, 0) Ki sin^ 9Qa^ + K2 sin^ 9q sin^ + 2X2 sin 9o cos 9q sin^ ^a. (9) 

The ground state of the FM particle with biaxial symmetry corresponds to the magnetization 
vector pointing in one of the two degenerate easy directions: 9 = 9o, and = 0,7r, other 
energy minima repeat the two states with period 27r. 

Performing the Gaussian integration over a, we can map the spin system onto a particle 
moving problem in one- dimensional potential well. Now the transition amplitude becomes 

/Ce = exp j-ii^ 1 - ^1 ^A^ cos6'o {(j)f - 4>i)^ J d4>exp{-SE[4>]) , (10) 



with the effective Euclidean action 



Se [(f>] = J 



(11) 



where X — m — -^y-, and U {(f)) — sin^ 9o sin^ cj). The potential U {(f)) is periodic with 
period tt, and there are two minima in the entire region 27r. We may regard the potential 
U {(f)) as a superlattice with lattice constant tt and total length 2?!, and we can derive the 
energy spectrum by applying the Bloch theorem and the tight-binding approximation. The 
translation symmetry is ensured by the possibility of successive 27r extension. 

Now we apply the periodic instanton method to evaluate the tunnel splittings of excited 
levels. The periodic instanton configuration (f)p which minimizes the Euclidean action of Eq. 
(11) satisfies the equation of motion 

where > is a constant of integration, which can be viewed as the classical energy of the 
pseudoparticle configuration. Then we obtain the kink-solution as 



= (12) 



sin (f)p — 1 — k sn {ljit, k) . (13) 



sn {cuiT, k) is the Jacobian elliptic sine function of modulus k, where 



and 



K2V sm Oq 



u;i^2^^KiK2smeo. (14b) 

In the case of resonant quantum tunneling at ground state with zero magnetic field, i.e., 
— > 0, A; — > 1, {u, 1) — > tanhu, A ^ 0, we have 

cos(f)p — tanh {ujit) (15) 

which is exactly the vacuum instanton solution derived in Ref. 21. 

The classical action or the WKB exponent can be obtained by integrating the Euchdean 
action Eq. (11) with the above periodic instanton solution. The result is found to be 



Sp = I dr 



+ 2EP, (16) 



with 



W = 2^/XS sin 00 [e {k) - (l - k'^) K (k) 



(17) 



where K {k) and E {k) are the complete elhptic integral of the first and second kind, re- 
spectively. In the low energy limit where E is much less than the barrier height, i.e., 
/c'^ = 1 - /c^ = HE / K2V sivl^ 9o < 1, we can expand K (k) and E (k) in Eq. (17) as powers 
of k' to include terms like k'"^ and k'"^ In {4:/k'), 

1 r. /4\ 1- 



E{k) 
K{k) 



1 + 



In 



In 



k' 



k'^ + 





1 










In 






+ 4 





k'^ + 



(18) 



With the help of small oscillator approximation for energy near the bottom of the potential 
well, E = E^'' = (n + I) Eq. (17) is expanded as 

1~ 



W = 2VXSsm9o - {n+-] + \ n 



n 



(19) 



23VA5sin^o 

Then the general formula Eq. (4) gives the low-lying energy shift of the nth excited level for 
FM particles with biaxial symmetry in a magnetic field along the hard anisotropy axis as 



23/2 y 



" y/KJll S 



^jKiK2sm9o (8v^,Ssin^o)"^'^'exp (^-2VXSsm9o) . 



(20) 



It is noted that Tt/S.S^'^ is only the level shift induced by tunneling between degenerate 
excited states through a single barrier. The periodic potential U (0) can be regarded as a 
one-dimensional superlattice. And the general translation symmetry results in the energy 
band structure, and the energy level spectrum can be determined by the Bloch theorem. It 
is easy to show that if f^'" are the degenerate eigenvalues of the system with infinitely high 
barrier, the low-lying energy level spectrum is given by the following formula with the help 
of the tight-binding approximation 



= - 2^A£r cos [tt (/. + 0] , 



?bia 



bia 



(21) 



where n = S 1 — (l + ^A^ cos^^o 1 and ^ is the Bloch vector which can be or 1 in the first 
Brillouin zone. Equation (21) includes the contribution of Wess-Zumino-Berry phase for 
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FM particles with biaxial symmetry at finite magnetic field. In the absence of a magnetic 
field, the tunnel sphtting is suppressed to zero for half-integer total spins by the destructive 
interfering Wess-Zumino-Berry phase. This topological quenching effect is in good agreement 
with the Kramers' theorem since the system has time-reversal invariance at zero field. In the 
presence of even weak external magnetic field this strict "selection rule" is relaxed, which 
leads to a finite tunnel splitting for half-integer total spins. The low-lying energy level 
spectrum is US^^" - 2hA£^''' cos{nS cosOq), and + 2^Ai'^" cos(7r^ cos ^o) for integer 

total spins. However, the low-lying energy level spectrum is hSn"" — 2hASn"' sin^irS cos 9o) , 
and hS^^^"" + 2hA£^"- sin(7rS' cos^o) for half-integer total spins. Therefore the tunnel splitting 
is /S.£n = 4A£^^*" |cos (7rS'cos6'o)| for integer spins, while the tunnel splitting is A£^„ = 
^/\gb%a I sin (vrS' cos 6*0) I for half-integer spins. The tunnel splitting will not be suppressed to 
zero even if the total spin is a half-integer at finite magnetic field. In Fig. 1, we plot the 
tunnel splitting in the magnetic field at the first excited level (n = 1) for integer total spin 
S = 100 by the analytical calculation and the exact diagonalization calculation, respectively. 
Here we take the typical values of parameters for single-domain FM nanoparticles Ki = 
lO^erg/cm^, A = 0.02, and the radium of particle r = 5nm. The analytical result is found 
to be in good agreement with the numerical result, which confirms the theoretical analysis. 
In Fig. 2, we plot the dependence of the first-excited-level splitting on the magnetic field 
for integer and half-integer total spins respectively, where the oscillation with the field and 
spin-parity effects are clearly shown. And in Fig. 3, we show the tunnel splittings of the 
ground-state level and the first excited level as a function of the magnetic field for integer 
total spins S = 100. It is clearly shown that the splitting is enhanced by quantum tunneling 
at the excited levels. 

Recently, spin systems have aroused considerable interest with the discovery that they 
provide examples which exhibit first- or second-order transition between the classical and 
quantum behavior of the escape rate.ii^ In general transitions in a metastable system can 
occur via quantum tunneling through the barrier and the classical thermal activation. It 
was showed that for a particle with mass m moving in a double- well potential U (x), the 
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behavior of the energy-dependent period of oscillations P (E) in the Euclidean potential 
—U {x) determines the order of the quantum-classical transition. If P {E) monotonically 
increases with the amplitude of oscillations, i.e., with decreasing energy E, the transition is 
of second order. The crossover temperature for the second-order phase transition is Tq^^ = 



ujo/2tt, ujq = \J\U" {xsad) \ I'm, wherc Xgad corresponds to the top (the saddle point) of the 
barrier, and ujq is the frequency of small oscillations near the bottom of the inverted potential 
—U [x). If, however, the dependence of P [E] is non-monotonic, the first order crossover 
takes place. 

Now we discuss the phase transition from classical to quantum behavior in this model. 
For the present case, the period of the periodic instanton is found to be P (E) = ^K{k). 
The monotonically decreasing behavior of P (E) is shown in Fig. 4, which yields that the 
second-order phase transition takes place. We found that the frequency of small oscilla- 
tions near the bottom of the inverted potential is uJi = 2^\/KiK2 sin Oq. The crossover 
temperature characterizing the quantum-classical transition is 



fcsTo^'^ = ^V^^i^2sin^o. (22) 

IT O 

At the end of this section, we discuss the possible relevance to the experimental test for 
spin-parity effects in single-domain FM nanoparticles. First we discuss the thermodynamic 
behavior of this system at very low temperature T ~ To = HAS^'^'^/kB- For FM particles with 
biaxial crystal symmetry at such a low temperature, the partition function of the ground 
state is found to be 

Z = 2 exp (-pnSi;''') cosh hf^fiAS^''' cos {ttS cos Oo)] , (23a) 



V 



for integer spins, while 

Z = 2 exp (^-(3hSl;''') cosh [2(3hASl;''' sin {nS cos ^o)] , (23b) 

for half-integer spins, where £^q*" = uJi/2. Then the specific heat is c = — T (9^JF/9T^), 
with JF = —kBTlnZ. For FM particles with biaxial crystal symmetry in the presence of a 
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magnetic field along the hard axis at sufficiently low temperatures, we obtain the specific 
heat as 

{cosh [2/3^A£o^^« cos (7r5 cos ^o)J I 

for integer spins, while 

c = 4ks ((3hAS-'^Y [sm{nScos0o)f _ ^24b) 

{cosh [2/?^A£:o^^" cos (tt-S cos ^o)] } 
for half-integer spins. In Fig. 5, we show the temperature dependence of specific heat at 

H/ He — 0.2 for integer and half- integer total spins respectively. It is clearly shown that the 

specific heat for integer spins is much different from that for half-integer spins at sufficiently 

low temperatures. 

When the temperature is higher tiASg"' <^ /c^T < hcui, the excited energy levels may 
give contribution to the partition function. Now the partition function is found to be 

, 2 



2o 



(25) 



1 + (l - e-'^^'^i) (y2l3hA8^'' cos (tt/x)) /q i^q^e-^^^'^ 

for both integer and half-integer spins. Zq = 2e~^^^l'^ j (^1 — e~^'^^^ is the partition func- 
tion in the well calculated for ksT <C AU over the low-lying oscillator like states with 
S^'' = (n + l/2)a;i. Io{x) = J2n=o {^f^f" / (r^-f is the modified Bessel function, and 
qi = 2^\/XS sin > 1- We define a characteristic temperature T that is solution of equation 
q^e~^^^'^^'^^ = 1. The temperature T = hu)i/2hiqi characterizes the crossover from ther- 
mally assisted tunneling to the ground-state tunneling. Then we obtain the specific heat up 
to the order of (phAEl^'^'f as 

c = kB {Phwif ^^^^ _ 2 + kB {y/2l3hA£^'' cos (tt/x)) {[2(1- e-^"^') + 4 {phwi) e'^"^- 
(e 1) 

- {iihw^f e-^^-^o (2gie-^'^"i/2jj _ ^^^^^^ [1 ^^-zfin^^ji _ 

x/^(2gie-^'^-^/2)}, (26) 

for both integer and half- integer spins, where Jg = — /i, and Jg = I2 — h/x. Iv{x) = 
Y,n=o (~1)" /n\r {n + u + 1), where F is Gamma function. 
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IV. MQC FOR TRIGONAL, TETRAGONAL AND HEXAGONAL SYMMETRIES 



In this section we will apply the method given in Sec. Ill to study resonant quantum 
tunneling of magnetization in single-domain FM nanoparticles with trigonal, tetragonal and 
hexagonal symmetry. For the trigonal symmetry, the anisotropy energy is 

E {e, (j)) = cos^ e-K2 sin^ 9 cos (30) - M^H cos 9 + E^, (27) 

where K2> 0. The energy minima of this system are at ^ = and = 0, |7r, |7r, and 

other energy minima repeat the three states with period 27r. This problem can be mapped 
onto a problem of one-dimensional motion by integrating out the fluctuations of 9 about 9o, 
and then the effective potential is 

t/(0) = ^sin3^osin2(^0). (28) 

Now U {(j)) is periodic with period |7r, and there are three minima in the entire region 27r. 
The periodic instanton configuration with an energy > is found to be 

sin^ (^'^p) = 1 - ^^^^^ (^2T, k) , (29) 

where k — y^'l — and U2 — ?>^/2^\/ K1K2 (sin ^o)^''^- The corresponding classical 

action is Sp — W + 2E(3, with 

95/2 

W = -^^/XS (sin 9of^^ [e (k) - (1 - k^) K {k)] . (30) 



The low-lying energy shift of the nth excited level is found to be 

6 V 
\pKn\ S 



n+1/2 

hASj;-' = -^-JKiK2 (sin 9of^ -^v^S (sin 9of^ exp —^/XS (sin 9of^ 



(31) 



The periodic potential U {(f)) can be viewed as a superlattice with lattice constant |7r and 
total length 27r. Then the Bloch theorem gives the energy level spectrum of the nth excited 
level S^' = (n+ 1^002 as 
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Ei" = h£t" - 2hA£t" cos 



(32) 



where $, = —1, 0, 1 in the first Brillouin-zone. The crossover temperature for the second-order 
phase transition is 



(33) 



For the tetragonal symmetry, 



E {6, 0) = cos' e + K2 sin^ 6 - K'^ sin* 6 cos (40) - MqH cos 6 + Eo, 



(34) 



where Ki ^ K2, K2 > 0. The energy minima are at 6* = 6*0 and = O,|7r,7r,|7r, and other 
energy minima repeat the four states with period 27r. The problem can be mapped onto a 
problem of particle moving in one-dimensional potential 



t/(0) = ^^sin^^osin^ {2(f)) 

f t 



(35) 



Now U (0) is periodic with period |7r, and there are four minima in the entire region 27r. In 
this case, the periodic instanton solution is 



sin (20p) = 1 — k sn {lv^t, k) , 



(36) 



where k 



HE 



2K'Vs\'a^ 6>o ' 



, and uj'i = A\/2^\/ K1K2 sin^ Oq. The associated classical action 



IS Sp^W + 2E/3, with 



W = 2^/^^/XS sin^ 00 [e (k) - (l - k') K (k) 



(37) 



The low-lying energy shift of the n-th excited level is 

8 V I - , ... ^ . x"+l/2 



" S 



^JKiK2 sin^ 9o (2^/'VXS sin^ Oq) exp (-2^/^VXS sin^ ^0) . (38) 



The periodic potential U (0) can be viewed as a superlattice with lattice constant ^tt and 
total length 2n. For this case the energy level spectrum of the nth excited level £^ — 

(n + Us is 
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Ei^ ^ hSi' - 2hASi' COS 



(39) 



where ^ = —1, 0, 1, 2 in the first Brillouin-zone. The crossover temperature is 

V 



For the hexagonal symmetry, the magnetocrystaUine anisotropy energy is 



(40) 



E {e, 0) = Ki cos' e + K2 sin^ 9 + sin^ 9 - K'^ sin^ 9 cos (60) - MqH cos 9 + Eo, (41) 

where Ki :$> K2, K^ K'^ > 0. The energy minima are at ^ = ^0 and (f) — 0, ^n, |7r, tt, |7r, |7r, 
and other energy minima repeat the six states with period 27r. Now the one-dimensional 
effective potential is 



[/(0) = ?^sin6^osin'(3</.). 



(42) 



For the present case, U (0) is periodic with period |7r, and there are six minima in the entire 
region 27r. The periodic instanton configuration is found to be 



sin (30p) — 1 — k sn {UJ4T, k) , 



(43) 



where k — ^1 — ^K'. v^in^e ~ 6-\/2^-\/^i^sin^ ^o- Then the corresponding classical 

action is obtained as Sp = W + 2£'/3, with 



93/2 

W = — v^^sin^^o 



E{k) - {l-k^)K{k) 



(44) 



Therefore the low-lying energy shift of the nth excited level is found to be 



3 X 22\/ 



'27/2 



n+1/2 



23/2 



/sTi /sTs sin^ ^0 VA5 sin^ ^0 exp — \f\S sin^ ^0 • (45) 



The periodic potential U (0) can be regarded as a one-dimensional superlattice with lattice 
constant ^tt and total length 27r. By applying the Bloch theorem and tight-binding ap- 
proximation, we obtain the energy level spectrum of the nth excited level = ('^ + l) ^4 

as 
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cos 



(46) 



In this case the crossover temperature characterizing the quantum-classical transition is 



In brief, the low-lying energy level spectrum of the magnetic tunneling states for trig- 
onal, tetragonal and hexagonal symmetry are found to depend on the parity of the total 
spins of the single-domain FM nanoparticles, resulting from the Wess-Zumino-Berry phase 
interference between topologically distinct tunneling paths. And the crossover temperature 
characterizing the quantum-classical transition is also obtained for each case. 



In summary, we have investigated the spin-phase interference effects in resonant quan- 
tum tunneling of the magnetization vector between excited levels for single-domain FM 
nanoparticles in the presence of a magnetic field along the hard anisotropy axis. The system 
considered in this paper has a general structure of magnetocrystalline anisotropy such as 
biaxial, trigonal, tetragonal and hexagonal symmetry. The low-lying tunnel splittings be- 
tween the nth degenerate excited levels of neighboring wells are evaluated with the help of 
the periodic instanton method in the spin-coherent-state path-integral representation. The 
low-lying energy level spectrum of the nth excited level is obtained by applying the Bloch 
theorem and the tight-binding approximation in one-dimensional periodic potential. This 
is the first complete study, to our knowledge, of spin-phase interference between excited 
levels and effects induced by magnetic field in FM particles with a general structure of 
magnetocrystalline anisotropy. 

One important conclusion is that for all the four kinds of crystal symmetries, the low- 
lying energy level spectrum depends on the spin parity significantly, resulting from the 
Wess-Zumino-Berry phase interference between topologically distinct tunneling paths. The 




(47) 



V. CONCLUSION 
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structure of low-lying tunneling level spectrum for trigonal, tetragonal or hexagonal symme- 
try is found to be much more complex than that for biaxial symmetry. The low-lying energy 
level spectrum can be nonzero even if the total spin is a half-integer for the trigonal, tetrag- 
onal, or hexagonal symmetry at zero magnetic field. Our study provides a generalization of 
the Kramers degeneracy to coherently spin tunneling at ground states as well as low-lying 
excited states. External magnetic field yields an additional contribution to the Berry phase, 
resulting in oscillating field dependence of the tunnel sphttings for both the integer and half- 
integer total spins. This oscillation effect can be tested with the use of existing experimental 
techniques. Due to the topological nature of the Berry phase, these spin-parity effects are 
independent of details such as the magnitude of total spins, the shape of the soliton and the 
tunneling potential. And the tunneling effect is enhanced by considering tunneling at the 
level of excited states. The transition from quantum tunneling to thermal activation is also 
studied. By calculating the oscillation period, we find the monotonically decreasing behavior 
of the period with increasing energy, which yields the second-order phase transition. The 
crossover temperature characterizing the quantum-classical transition is obtained for each 
case. Comparison with the numerical results gives strong support for the analytical calcu- 
lations presented in this paper. The heat capacity of low-lying magnetic tunneling states is 
evaluated and is found to depend significantly on the parity of total spins for FM particles at 
sufficiently low temperature. This provides a possible experimental method to examine the 
theoretical results on spin-phase interference effects. Our results presented here should be 
useful for a quantitative understanding on the topological phase interference or spin-parity 
effects in resonant quantum tunneling of magnetization in single-domain FM nanoparticles. 

The theoretical calculations performed in this paper can be extended to the single-domain 
antiferromagnetic nanoparticles, where the relevant quantity is the excess spin due to the 
small noncompensation of two sublattices. Work along this line is still in progress. We hope 
that the theoretical results obtained in the present work will stimulate more experiments 
whose aim is observing the topological phase interference or spin-parity effects in resonant 
quantum tunneling of magnetization in nanoscale single-domain ferromagnets. 
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Fig. 1. The tunnel splitting Asi for biaxial FM particles at the first excited level 
(n = 1) as a function of the field H for integer spins {S — 100) by the analytical and the 
exact diagonahzation calculation, respectively. Here Ki — lO^erg/cm^, A = 0.02, and the 
radium of particle r — 5nm. 

Fig. 2. The tunnel sphtting Aei for biaxial FM particles at the first excited level {n — 1) 
as a function of the field H for integer {S = 100) and half-integer {S = 100.5) total spins 
respectively. Here Ki = lO^erg/cm^, A = 0.02, and the radium of particle r = 5nm. 

Fig. 3. The tunnel splittings A£„ for biaxial FM particles at the ground level {n — 0) 
and the first excited level (n = 1) as a function of the field H for integer {S = 100) total 
spins respectively. Here Ki = lO^erg/cm^, A = 0.02, and the radium of particle r = 5nm. 



Fig. 4. The relative period P (E) /P [E — K2V s\n^ Oojol the periodic instanton as a 
function of energy E/K2Vfim^ 00 in the domain < hE < K2V sin^ ^^0- 

Fig. 5. The low-temperature specific heat for biaxial FM particles with integer {S = 100) 
and half-integer {S = 100.5) total spins respectively. Here A = K2/K1 = 0.02, H/H^ = 0.2. 
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